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We calculate the local density of states of a two-dimensional electron system under strong crossed 
magnetic and electric fields. We assume a strong perpendicular magnetic field which, in the absence 
of in-plane electric fields and collision broadening effects, leads to Landau quantization and the well- 
known singular Landau density of states. Unidirectional in-plane electric fields lead to a broadening 
of the delta-function-singularities of the Landau density of states. This results in position-dependent 
peaks of finite height and width, which can be expressed in terms of the energy eigenfunctions. These 
peaks become wider with increasing strength of the electric field and may eventually overlap, which 
indicates the onset of inter-Landau-level scattering, if electron-impurity scattering is considered. 
We present analytical results for two simple models and discuss their possible relevance for the 
breakdown of the integer quantized Hall effect. In addition, we consider a more realistic model for 
an incompressible stripe separating two compressible regions, in which nearly perfect screening pins 
adjacent Landau levels to the electrochemical potential. We also discuss the effect of an imposed 
current on the local density of states in the stripe region. 

PACS numbers: 73.43.Cd 
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I. INTRODUCTION 

The integer quantum Hall effecti (IQHE) is one of the 
most important discoveries of condensed matter physics, 
observed on two-dimensional electron systems (2DES) 
subjected to a strong magnetic field B perpendicular to 
the plane of the system. Measuring the low-temperature 
resistance in such systems, one finds that in certain B- 
intervals, the "plateau regimes" of the IQHE, the longitu- 
dinal resistance vanishes, indicating dissipationless trans- 
port, while the Hall resistance assumes quantized values, 
Rh = h/{ve 2 ), where h = 2nh is Planck's constant, 
e the elementary charge, and v a positive integer. An 
idealized homogeneous 2DES in a heterostructure (like 
GaAs/AlGaAs) without any scattering or interaction of 
the electrons, subjected to a homogeneous in-plane elec- 
tric field, would yield the same resistance values. For this 
ideal 2DES, v = 2nl 2 n c i is the filling factor of the Lan- 
dau levels, with energy eigenvalues E n = hw c (n + l/2) for 
n = 0, 1, . . . , n c \ is the electron density, I = y/ch/e\B\ 
the magnetic length (we use CGS units, with c the ve- 
locity of light), and uj c = e\B\/(m*c) the cyclotron fre- 
quency, with m* the effective mass (m* = 0.067m e , for 
GaAs). In this paper we neglect spin-splitting of the Lan- 
dau levels and take account of the spin degree of freedom 
by a degeneracy factor g s = 2. The fact, that in the 
IQHE one observes integer values of v, suggests that, in 
the plateau regime, the electrons which carry the cur- 
rent fully occupy an integer number of Landau levels, 
and thus are in states which are energetically separated 



from empty states, so that at low temperatures dissipa- 
tive scattering processes are suppressed. Then one has 
to understand, why the relevant filling factor v remains 
constant over a S-interval of finite width, and which pro- 
cesses limit the plateau regimes and lead back to dissi- 
pative transport. 

If the effective filling factor has to be constant over 
an interval of B-values, obviously the relevant electron 
density must change. In the early attempts to ex- 
plain this phenomenon, it has been assumed that the 
Coulomb interaction between the electrons is unimpor- 
tant for the understanding of the IQHE. Localization 
theories^, stated that electron-impurity interaction leads 
to localized states in the tails of collision-broadened Lan- 
dau levels, which do not contribute to the current trans- 
port. The existence of such inert localized states may ex- 
plain the resistance quantization, even if a homogeneous 
current distribution is assumed. Such theories have, how- 
ever, problems with explaining the enormous accuracy 
of better than 10 -8 , with which the quantized resistance 
values can be reproduced^ even in narrow Hall bars with 
a width of a few micrometers^ The confinement of the 
2DES to the interior of such Hall bars leads to edge states 
in the gaps between bulk Landau levels, which also have 
been considered to be important for the transport in the 
plateau regime of the IQHEi&i This picture, which leads 
to extremely high current densities in the edge states and 
strongly independent electron density profiles, has been 
critizised 8 - because it neglects important screening effects. 

The Landau quantization, which alters the energy- 
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independent density of states (DOS) of the 2DES at 
B = 0, D = m*/(nh 2 ), to the Landau DOS with sharp 
peaks around the Landau energies E n , leads to pecu- 
liar screening effects, with nearly perfect screening if the 
Fermi energy Ep (or the chemical potential) coincides 
with a Landau energy, and with no screening, if Ef falls 
into a gap between two Landau levelsAi^ In an inho- 
mogencous 2DES, as e.g. in a Hall bar with electron 
depletion near the edges, this should lead to the occur- 
rence of "compressible regions" , where a Landau level is 
pinned to the Fermi energy and the density is position- 
dependent, and "incompressible regions" which separate 
neighboring compressible regions with adjacent Landau 
levels at Ef- In these incompressible regions Ef lies 
in the gap between two Landau levels, and one expects 
there a constant electron density corresponding to an in- 
teger value of the filling factor.— itr— Scanning force mi- 
croscope experiments 1 ^ - — on narrow Hall bars have con- 
firmed this picture. In the plateau regimes of the IQHE, 
the Hall potential across the sample drops across stripes 
at the positions predicted for incompressible stripes (ISs), 
and is constant elsewhere. Well outside the plateau 
regimes, the Hall potential varies nearly linearly be- 
tween the sample edges. Self-consistent calculations for 
equilibriu m 12 ' 13 and transport 1 ^— have clarified the sit- 
uation. At sufficiently high temperature the conductivity 
is Drude-like, the current density j y (x) along the sample 
is proportional to the electron density n e \{x), and the 
Hall potential varies linearly across the sample. With 
decreasing temperature, near the lines of constant local 
filling factor v(x) = 2'Kl 2 n \(x) — kg s with even-inter 
values (since spin-degeneracy is assumed) local minima 
of the longitudinal resistivity develop (Shubnikov-deHaas 
effect), which lead to local maxima of the current density. 
With further decreasing temperature, ISs with extremely 
small values of the longitudinal resistivity develop along 
these lines and the current becomes confined to these 
stripes, where it can flow nearly without dissipation. It 
has been argue d 18 ' 19 that an IS supporting dissipationless 
transport between two compressible regions can develop 
only if the distance between these regions is sufficiently 
large, larger than several times the extent of a typical 
wavefunction, say about 7£, which requires a sufficiently 
slow variation of the effective potential^ For the consid- 
ered class of samples this guarantees that, in agreement 
with the experiments, the plateau regimes of the IQHE 
are well separated on the J3-axis, and that only ISs with 
the same value of the local filling factor can exist simul- 
taneously in the Hall bar. 

If the stripe between two compressible regions becomes 
very small, one may expect that energy eigenfunctions 
with centers in different regions overlap and lead to quasi- 
elastic inter-Landau- level scattering (QUILLS), a mech- 
anism which has been discussed for a long time as a 
possible reason for the breakdown of the IQHE under 
strong imposed currentsi 21 ' 22 Such an overlap of wave- 
functions, belonging to different Landau quantum num- 
bers and having different center coordinates, but having 



the same energy eigenvalues, can occur in Landau lev- 
els tilted by a constant (Hall) electric field. A suitable 
quantity to study such overlap effects seems to be the 
local density of states (LDOS), which has recently been 
studied for such a constant-electric-field model and in- 
terpreted with respect to the IQHE ) 23 - 24 however with a 
rather complicated mathematical approach and dubious 
results for the IQHE. 

The purpose of this paper is to present a simple way 
to calculate the LDOS for a 2DES with translation sym- 
metry in one in-plane direction, sect|TTl to give explicit 
analytic and numeric results for simple model potentials 
in the other in-plane direction, sect JIIH and, finally, to 
evaluate the LDOS for a simple but realistic model of 
an IS between two compressible regions, which carries 
an intrinsic and, possibly, an imposed external current, 

sect us 

II. ELECTRIC-FIELD-BROADENED LANDAU 
LEVELS 

We describe a 2DES in the x-y-plane, subjected to a 
strong magnetic field B = (0,0,5) = V x A(r) in z- 
direction, in an effective-field (e.g. Hartree) approxima- 
tion by a single-particle Hamiltonian 

where the potential energy V(r) may contain the effect of 
externally applied static electric fields, of lateral confine- 
ment, and of the average Coulomb interaction with the 
other electrons of the 2DES. Once the eigen-functions 
ip a (r) of the Schrodinger equation 

(H-E a )Mr) =0 (2) 
are known, one can calculate the electron density 

n(r) = £/ Q |^ Q (r)| 2 , (3) 

a 

where the occupation probability f a of the energy eigen- 
state | a) may depend on all the quantum numbers of 
conserved quantities collected in a, i.e., two for orbital 
motion and one for spin. 

A. Local density of states (LDOS) 

If in Eq. Q the occupation probability of the state \a) 
depends only on its energy eigenvalue, f a = f(E a ), it 
may be useful to express the density 

n(r) = J dE f(E) D(E; r) (4) 

in terms of the "local density of states" (LDOS): 

D(E;r)=J2t{ E - E <*)\M*)\ 2 - (5) 

a 
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This formula for the LDOS is easily generalized to include 
the effect of quasi-elastic scattering of the electrons by 
randomly distributed impurities, which leads to a "col- 
lision broadening" of the (5-function in Eq. (|5|). System- 
atic calculations of collision broadening^— usually start 



from the Green operator Gi mp (z) = (z — H — V u 



for 



a fixed impurity configuration, described by an impu- 
rity potential Vi m p(r), and calculate approximately the 
average over all possible impurity configurations. This 
average is expressed in terms of a self-energy operator 
E(z), 

G(z) = [z-H -S(z)]- 1 = (G imp {z)) . (6) 

\ / imp 

With G + (E) = G(E + i0 + ) the corresponding general- 
ization of Eq. (J5J) reads 



D(E;r) = Im(r|G + (S)|r) 



(7) 



Without impurities £(z) = 0, and Eq. (|Tj) reduces to 
Eq. ([5]), with the notation (r|a) = i/j a (r). 



B. Translation symmetry in y-direction 

In the following we assume that the system is 
translation-invariant in y-direction, but electric fields in 
^-direction, E = (E x ,0,0) = VV(x)/e, will be allowed. 
The translation invariance in y-direction suggests the 
Landau gauge A(r) = (0,xB,0) for the vector poten- 
tial, so that the single-electron Hamiltonian ([1]) becomes 
cyclic in y and allows the separation ansatz 



ip(x,y) 



Jky 



(8) 



where L y (— > oo) is the normalization length in y di- 
rection, and the quasi-continuous momentum quantum 
number k assumes the values k = 2Trn y /L y , for arbitrary 
integers n y . With this ansatz the Schrodingcr equation 
([2]) reduces to the one-dimensional form 



Hx<p n ,x(x) = E n (X)cp n> x(x), 



(9) 



with the effective Hamiltonian 



H 



where X = —£ 2 k denotes the center of the parabolic po- 
tential, which describes the effect of the magnetic field 
and leads for fixed X to a discrete energy spectrum 
E n (X). Here and in the following we neglect spin split- 
ting and consider spin by a degeneracy factor g s = 2. 
In general the eigenstates (r\n,X) carry current in y- 
direction, and the expectation value of the velocity oper- 
ator v y is given by (Hcllmann-Feynman theorem) 



(n,X\v y \n,X) 



1 dE n {X) 



m^uJc 



dX 



1 dE r , 



h dk 



(11) 



Then the electron density, Eq. ((3|), depends only on x, 



n{x) 



9« 



2ir£ 2 



Y,JdXf n , x \ Vn , x (x)\ 2 , (12) 



and is accompanied by a current density 

dE n (X) 



3y 



dX 



ip n ,x(x)\ ■ (I 3 ) 



The LDOS, Eq. ©, reduces to 
D(E; X )=to^r dX6 {E-E n {X)) |^, x (x)| 2 . (14) 

If the dependence of E n (X) on X is smooth enough to 
allow for a Taylor expansion around the center coordi- 
nate X U: e defined by E n (X n ,E) = E, the X-integral in 
Eq. (TT4")) can be evaluated: 



D{E;x) 



9s Wn,x n , E {x)\ 2 
2ir£ 2 4- \E' n (X ntE )\ ' 



(15) 



with E' n (X nt E) = dE n /dX(X n: E)- Before we illustrate 
some properties of this LDOS with typical examples, we 
introduce a simple treatment of collision broadening. 



C. Collision broadening 

1. Homogeneous 2DES without electric field 

For V(x) = we get the well known Landau problem 
with energy eigenvalues and eigenfunctions 

E n = hbj c (n+ i), <p n ,x(x) = -^=u n ( X Y (16) 



respectively, where the normalized oscillator wavefunc- 
tions, 



1 \V2 



2 n nU/ir 



-C 2 /2 



(17) 



29 



are given by the Hermite polynomials H n (£) of order : 
Since here the energy eigenvalues are independent of X , 
the X-integral in Eq. (|T4l) reduces to the normalization 
integral of the eigenfunctions, and the LDOS reduces to 
the well known Landau DOS of the homogeneous system 



D(E;x) 



9s 



2tt£ 2 



(18) 



which does not depend on the position x. To include 
the effect of collision broadening, one has to evaluate 
the self-energy operator. With weak assumptions (like 
rotation symmetry) on the impurity potentials, one can 
show that E(z) and the Green operator G(z) are diagonal 
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in the Landau representation, and that the matrix ele- 
ments together with the eigen- energies E n (X) do not de- 
pend on XM~2& Then in Eq. (JTHJ> the singular 5(E — E n ) 
is replaced by a spectral function A n (E — E n ) of fi- 
nite width. Depending on the approximation scheme, 
several analytical forms for the spectral function have 
been obtained. The self-consistent Born approximation 
(SCBA) 28 ' 30 leads, if scattering between different Landau 
levels is neglegted, to a semi-elliptical form, 



A s n CBA (E-E n ) 



1 



tiT„ 



1 - 



E — E n 



2L„ 



while other approaches yield a Gaussian form 
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A%(E-E n ) 



1 



2?r r, 



exp 



E — E„ 



(19) 



(20) 



In the limit of short-range impurity potentials the ma- 
trix elements of the self-energy and thereby the L„ in 
Eqs. (fl~9| and (|20|) become even independent of the Lan- 
dau quantum number n. 



III. EXACTLY SOLVABLE MODELS 
A. Constant electric Held 

Simple analytic results are also obtained for the case of 
a constant in-plane electric field E = (£^,0,0), leading 
to the potential V(x) = exE x . Within classical mechan- 
ics, this leads for an ideal 2DES to a constant Hall drift of 
the centers of the cyclotron motion, which can be elim- 
inated by a Galilei transformation to a coordinate sys- 
tem moving with the drift velocity vrj = cE x B/B 2 = 
(0, — cE x /B, 0). Since all electrons suffer the same drift 
velocity, the current density = -evflnfi) is pro- 
portional to the electron density n(x), and one obtains 
Ohm's law j(x) = ct(.t)E with the Hall conductivity 
Cyx(x) = (ec/ B)n(x) and vanishing longitudinal conduc- 
tivity, a xx (x) = 0. 



1. Eigenstates and LDOS 



2. Model for non-homogeneous systems 

We now consider the more general case that the ef- 
fective Hamiltonian (fT0|) contains a position-dependent 
potential V(x). For each value X of the center coordi- 
nate this leads again to a discrete energy spectrum, but 
these generalized Landau energies usually depend on X 
and form dispersive Landau bands with energies E n (X). 

As mentioned above, in the absence of electric fields 
the Landau energy eigenvalues do not depend on the 
center coordinate X, and as a consequence Green and 
self-energy operator are diagonal in the Landau represen- 
tation and independent of X. In the presence of electric 
fields, however, the eigen-energies depend on X, and so 
do Green and self-energy operators. Moreover we can 
show that both are no longer diagonal in the generalized 
Landau representation, which diagonalizes the Hamilto- 
nian in the absence of collisions. This more complicated 
situation is found even in the case of short-range im- 
purity potentials, which for homogeneous systems with 
V(x) = leads to the same collision broadening for all 
Landau levels. Then the evaluation of collision broaden- 
ing effects becomes much more complicated, and exceeds 
the scope of the present paper. 

For weak electric field E x , on the other hand, the for- 
mula 

D(E-,x)=^^fdXA n (E-E n {X))yl :X {x), (21) 

n 

which is correct for E x = and a straighforward general- 
ization of Eq. ([5]), should yield a reasonable description of 
collision broadening effects. Since we are not aware of a 
better, practicable approach to treat collision broadening 
in the presence of electric fields, we will in the following 
use Eq. ([21]) as a phenomenological rule to estimate the 
consequences of such scattering effects. 



Inserting V(x) = exE x into the Hamiltonian (|10[) re- 
sults in a shifted parabolic potential with the new center 
X = X — eE x / (m* ui 2 ) and position-independent terms, 
which add to the oscillator energies e„ = hw c (n + 1/2). 
The resulting energy eigenvalues and eigenfunctions are 

E„{X) = £„ + eE x X + ^-v 2 D 

= e„ + eE x X - —v 2 D = E n (X), (22) 

and 

<Pn,x(x) = ^q u ™( X £ )> ( 23 ) 

respectively, with vd = cE x /B. From Eq. (fTTj) 
we see that each state carries the same current 
—e{n, X\v y \n, X) = e 2 E x /m*u! c = evo, in analogy to 
the fact, that the radius of the classical cyclotron orbit 
has no influence on the drift velocity of its center. As a 
consequence of Eqs. (TT2")) and (fT3"|) the current density is 
directly proportional to the electron density, 

j y {x) = ev D n(x), (24) 

independent of the occupation probability of the eigen- 
states, just as in the classical case. 

Due to the linear dependence of E n (X) on X, Eq. (fTS"]) 
can be written as 

nfj? ^ 9s \ - 1 2 ( E n {x) - E ^ 
g( ^* ) = 2^4^M eE x i )■ (25) 

This result has been obtained in Ref. [23| in a much less 
transparent way, starting from the symmetric instead of 
the Landau gauge for the vector potential. 

Since in Eq. (|25[) the argument of the wavefunctions 
depends linearly on both, the position x and the energy 
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FIG. 1: (color online) Heavy black lines (7 = 0): LDOS 
according to Eq. (|25p for three values of the electric field 
strength, \eE x \£/ftoj c = 0.1, 0.2 and 0.3. Also shown is the 
effect of level broadening according to Eqs. ()2ip and (|20p for 
7 = T/hw c = 0.05 and 0.1. D = m'/^h 2 ), g s = 2. 



E, the energy dependence of the LDOS reflects the posi- 
tion dependence of the energy eigenfunctions. The con- 
tributions of the individual Landau levels, which become 
(5-function like for vanishing E x , become wider with in- 
creasing electric field, with a width proportional to eE x £. 
If one shifts the position from itoi + o, one obtains the 
same profile for the LDOS, but shifted on the energy 
axes by —eE x a: D(E; x + a) = D{E — eE x a; x). Typical 
results 2 ^ for the LDOS according to Eq. (f25j) are shown 
by the heavy black lines in Fig. [TJ It is seen that the 
gaps between the lowest adjacent Landau levels close for 
0.1 < eE x £/hu> c < 0.2, whereas the zeroes of the LDOS 
determined by the zeroes of the eigenfunctions remain 
rather stable. 

If contributions to the LDOS due to adjacent Landau 
levels overlap, this means that, at the same position and 
at the same energy wavefunctions due to different Lan- 
dau levels have finite values. Then, if there is also a 
non- vanishing impurity potential at this position, the po- 
tential matrix element between these adjacent levels is 
finite and there must be elastic scattering between these 
levels. This quasi-elastic inter-Landau-level scattering 
(QUILLS) has been discussed for a long tim o 21 ' 22 as a 
possible mechanism for the breakdown of the IQHE un- 
der strong imposed currents. Apparently QUILLS must 
become important in the neighborhood of narrow incom- 
pressible strips, where the local potential must bridge an 
amount of order tvjj c across a strip of a width less than 
about 101 

The colored lines in Fig. [T] are calculated from Eq. (f2"Tj) 
for a Gaussian spectral function, Eq. (|2"0"|) . with n- 



independent r„ = 7 hio c for two values of 7. Appar- 
ently the zeroes of the LDOS, which are due to the ze- 
roes of the energy eigenfunctions, are smeared out al- 
ready by a very weak collision broadening, and are of no 
importance in real samples. A discussion 2 ^ of a possi- 
ble importance of these zeroes for the QHE is therefore 
without any relevance. On the other hand, the value of 
the LDOS in the gap between two adjacent Landau lev- 
els is of importance. In order to yield a plateau in the 
IQHE, the gap in an incompressible strip between two 
adjacent compressible regions must be sufficiently well 
developed. As a measure for the quality of such gaps 
we may consider the overlap of the contributions of adja- 
cent Landau levels to the LDOS, according to Eq. (|25p . 
We define the overlap as the product of these contribu- 
tions in the middle E n _ n+ i(x) = [E n (x) + E n +i(x)]/2 be- 
tween these levels, devided by the square of the zevo-B 
DOS Do = m*/(Trh 2 ), to make the overlap dimension- 
less. Since g s /(2Tr£ 2 D ) = hu c , Eqs. ^ and yield 
for the dimensionless overlap of level n and n + 1 : 

0„,„+i(r/) ^ -^u 2 (- ^-)w 2 +1 (^-), (26) 

with -q = e\E x \ijhbj c . The results for the lowest 
gaps, Oo,i(»?) = exp(-l/2?7 2 )/(27r?7 4 ) and Oi )2 (r?) = 
O ,i(f?)(2 - l/V) 2 /8, are plotted in Fig. HJ If we say 




il 



FIG. 2: Dimensionless overlap for the two lowest Landau gaps 
as function of f] = e\E x \l/ftw c . 

that the gap between Landau level n and n + 1 is well 
developed if O n , n +i < 10~ 8 , this defines a critical value 
Vnn+i (Vo~i ~ 0-15, rif 2 ~ 0.13) and thereby a critical 
field-strength E^ r n+1 = rj^ n+1 tiij c /e£. Only for suffi- 
ciently small electric fields with 

\E X \ < 2.1 ?7 - 1+1 (B/10T) 3 / 2 x 10 6 V/m (27) 

the gap between the Landau levels n and n + 1 is well de- 
veloped. An equivalent formula, with 2.1?7^ r n+1 replaced 
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by l/[V2n + l+V2n + 3j, was given in Eq. (41) of Ref.l23l 
as experimentally verified criterion for the breakdown of 
the IQHE. 

2. Occupation of eigenstates 

The properties of a 2DES arc not only determined by 
the properties of the single-particle eigenstates, but also 
by their occupation. Here we consider two different ex- 
amples. 

According to the general rules of statistical mechan- 
ics, a thermal equilibrium state is characterized by the 
expectation values of its conserved quantities, which in 
our case are (1) the particle number, or for fixed vol- 
ume the average particle density, (2) the energy, and (3) 
the quasi-momentum in y-direction, p y with eigenvalues 
hk = —m*Lo c X, related to the translation symmetry. The 
grand canonical distribution function under these bound- 
ary conditions yields 

fn.x = fW[E n (X) — fi— hv y X/e 2 }), (28) 

with f(x) = 1/(1 + e x ) the Fermi-Dirac function and 
8 = 1/fcsT, — and —(3hv y /£ 2 Lagrange multipliers 
conjugated to the conserved quantities energy, particle 
number, and quasi-momentum, respectively. 

a. Constant electron density. If we choose v y = vd, 
equal to the classical Hall drift velocity, the argument of 
the distribution function becomes independent of X, and 
the JT-integral in Eq. (fT2j) reduces to the normalization 
integral of the wavefunctions, so that the electron density 
becomes independent of the position x and is given by 

n(x) ^d = ^£ /GSfc, - m*v%/2 - M ]), (29) 

n 

which apart from an unimportant shift of the energy zero 
is the same as in the absence of the electric field E x . This 
choice of the Lagrangian multiplier apparently describes 
the quantum analog of the homogeneous Hall system, 
with spatially constant electron and current densities. 

Considered as a function of the chemical potential, 
the electron density n c i(/i) for fixed magnetic field is a 
step-function, with steps of height g s /(2Tr£ 2 ) and a width 
of the order of fc^T (or of T n if collision broadening is 
considered), located near the Landau energies fj, ~ e n . 
At low temperatures, ksT <C huj c , (and for weak colli- 
sion broadening, T n <C hu> c ) these steps are separated 
by wide plateaus of constant n c i(/i), where [i varies in 
the gap between two adjacent Landau levels (LLs). In 
these plateaus there are no states at the Fermi energy 
(i.e. near /z) and no quasi-elastic scattering is possible. 
Therefore, the longitudinal conductivity is zero and the 
Hall conductivity has the quantized value an = e 2 v/h 
with an integer value of the filling factor v = 27r^ 2 n G i(/i), 
even if one allows for quasi-elastic impurity scattering, 
which may lead to dissipation if /z is located in a broad- 
ened LL. This scenario is considered in Ref. [23| in order 



to explain the IQHE. The same physical situation can 
be considered for constant /i and varying magnetic field 
B. With increasing B both the degeneracy g s /{2 / n£ 2 ) of 
the LLs and their energies e n increase linearly with B. 
If fi is within a temperature and collision broadened LL, 
we may observe dissipation and the conductivity compo- 
nents are not quantized. With increasing B the energy 
of this broadened level rises above the chemical poten- 
tial, which then falls into the gap below this broadened 
level. Then the electron density increases at constant 
filling factor linearly with B, and the conductivity com- 
ponents are quantized, until the next lower broadened LL 
reaches the energy [i. The dissipation sets in again, and 
the filling factor decreases by g s as this level passes the 
fixed chemical potential, which causes a rapid decrease 
of the electron density. 

If in the low-temperature transport experiments on 
2DESs the chemical potential would be constant, this 
scenario would explain the IQHE. This scenario, which 
at low filling factors requires that changing the magnetic 
field induces large density changes (up to 50%), is, how- 
ever, unrealistic. In real experiments such a large elec- 
tron exchange between the 2DES and its surrounding 
is hardly possible, and the assumption of constant elec- 
tron density is much more realistic than that of constant 
chemical potential. For constant n e i(fi), Eq. (|29[) requires 
that /I oscillates as a function of B, and these oscilla- 
tions have indeed been observed experimentall y 32 ' 33 If a 
broadened LL is completely occupied, a further lowering 
of B at constant ri e i(/x) requires, that /x jumps to the 
lower edge of the next higher LL, so that the quantized 
values of the conductivity components occur only at a 
single, isolated value of B, not in a whole B-interval. 
This explains the well-known Shubnikov-deHaas effect, 
but not the IQHE as is claimed in Ref. [23j |. The at- 
tempts of Refs. HU and [24[ to calculate the conductivity 
tensor are also not compatible with accepted transport 
theorie a 28 ' 30 and yield incorrect results. 

If one writes in Eq. (J2"g| /i*(A) = fj, — hv y X/£ 2 and re- 
places the center coordinate by the position x, one gets a 
position-dependent electrochemical potential n*(x). This 
may reasonably describe a stationary non-equilibrium, 
and possibly dissipative, state as is found in a quantum 
Hall system outside the plateau regime of the QHE i 17 i 18 
In a thermal equilibrium state, however, the electrochem- 
ical potential /i* must be spatially constant and can de- 
pend only on eigenvalues of conserved quantities, such as 
X and E n {X). 

A 2DES with homogeneous electron and Hall current 
densities may be a good approximation to the interior of a 
laterally confined system, far away from the edges. The 
same single particle states and energies may, however, 
also be used to describe an edge region of a laterally 
confined system. 

b. Variable electron density. Let us now focus on an 
edge region and assume that the confinement potential 
there can be approximated by the linear potential V(x) = 
eE x x. Let us further assume that the total, laterally 
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confined, system is in thermodynamic equilibrium, with 
vanishing total current. To describe such a state, we put 
in Eq. (|28|) v y — and take a fixed constant value for 
the electrochemical potential. Then the electron density 
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FIG. 3: (color online) Upper part: filling factor v(x) — 
2tv£ 2 ti(x) for fi* = 1.5huj c , ksT = 0.0lhu c ; lower part: poten- 
tial V(x) = eE x x, energies E n (X) (dashed lines) and the 
wavefunctions squared at X 3 — and X s = ±fkJ c /eE x , 
shifted by their eigenenergies in units of the cyclotron en- 
ergy i? C yc = ftu) c . In both parts 77 = eE x l/ftw = — 0.1. The 
dotted lines indicate the spatial extent of the wavefunctions. 

can be calculated directly from Eq. ()12l) or, cquivalcntly, 
from the LDOS (2SD, 

^)^E/f/(^^K.(^) 

= j dEf[^-f)D(E;x). (30) 

In the lower part of Fig. |3] we show, for E x < 0, the 
linear potential V(x) = eE x x, the lowest energy eigenval- 
ues E n (X), and, for X = the squared eigenfunctions 
ul(x,X) = ul([x - X]/£), as defined in Eq. flT7|), but 
shifted upwards by the amount of their energy eigen- 
value. Also shown are the squared eigenfunctions at 
X± = ±fno c /eE x , where E n (X±) = E n ±i(0). Obvi- 
ously the distance \X±\ between the center coordinates 
of eigenfunctions, which have the same energy but be- 
long to adjacent Landau levels, is inversely proportional 
to the field strength \E X \, whereas the shape of the wave- 
functions is independent of E x . Their spatial extent, 
\x — X\ < 1.2R„, is indicated by the dotted lines in Fig. 02 
which we estimate by the Landau radius R n = t\/2n + 1 
defined by (m*/2) uj 2 R 2 n = hui c (n + 1/2). Since the dis- 
tance between the eigenfunctions shrinks with increas- 
ing \E X \, their overlap increases, and is larger for higher 
than for lower Landau quantum numbers. This increas- 
ing spatial overlap of wavefunctions with the same energy 



eigenvalue has, of course, the same origin as the increas- 
ing energetic overlap in the LDOS D(E; x) at a fixed 
position x. 

The upper part of Fig. [3] shows the electron density 
calculated from Eq. (|30| for /i* = 1.5hui c , i.e., the energy 
for which the eigenfunctions with n = 0, 1, and 2 are in- 
dicated. Apparently n(x) increases with a; in a stepwise 
manner, with steps at positions where wavefunctions of 
states with energy E n (X) < /1* become relevant. The 
extent of the wavefunctions determines the width of the 
steps, and their zeroes lead in the limit T — > to zeroes 
of the slope dn/dx. This internal structure of the steps 
occurs on a length of the order of the magnetic length 
I (~ lOnm for B ~ 10T), and has never been resolved 
in real Hall bars (with a width > 10 /im). The width 
of the plateaus between the steps increases inversely pro- 
portional to \E X \. We want to emphasize that a change of 
the electrochemical potential /i* affects the density profile 
only by a rigid shift, n(x; (i* + <5/i) = n(x — S(i/eE x ; [i*). 

In the limit T — > the density is given by 

n(x;E F )= dE D(E; x), (31) 
J — 00 

with Ep — /i*(T = 0). On the scale of Fig. [3] this cannot 
be distinguished from the given result for ksT = 0.0lhuj c . 

Of course Eq. (|24|) yields for the Hall conductivity the 
trivial result an (x) = (ec/B)n(x; Ep), in agreement with 
Eq. (22) of Ref. [24[ (which considers j x and E y instead of 
our j y and E x ). To model the current through a macro- 
scopic device by I y = L dxj y {x) with j y (x) = uh{x)E x 
calculated from the present constant- E x model, as is done 
in Ref. [24[, is not meaningful since it effectively intro- 
duces a very unphysical description of the sample edge 
at x = W . Describing this edge by a reasonable con- 
finement potential, and the state of the system by the 
Fermi-Dirac distribution with a constant Fermi energy, 
would lead to vanishing total current, I y = 0. To de- 
scribe a dissipative state with non-vanishing total cur- 
rent, one needs a position-dependent electrochemical po- 
tential. A dissipation-free state with finite total current 
can be described as we explained above, but not with a 
position-independent distribution function that depends 
only on energy. 

In Ref. [23[ n(x; Ep) was calculated and discussed for 
x = as function of Ep. It was speculated that the 
structure of this curve might be related to the quantized 
Hall effect. This is, however, incorrect. The QHE is 
observed on real, confined systems, where a finite total 
current, and voltages along and across the sample, can be 
measured. The resistance quantization is a property of 
the sample as a whole. It can not be explained by local 
properties like energy dependence of the electron den- 
sity at a position x somewhere inside the sample. The 
present linear-potential model, on the other hand, leads, 
if taken serious, to an electron density n(x; Ep) and a cur- 
rent density j y (x) = (ecE x j B)n(x; Ep), which increase 
with increasing x, just because the number of eigenstates 
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with energy eigenvalues E n (X) < Ep increases with in- 
creasing X. This model makes sense only as an approxi- 
mation to an edge region of a laterally confined sample, 
which then, with the choice of a constant electrochemical 
potential, has vanishing total current. The property of 
the n{x = 0; Ep)-cmve tells nothing about such a real 
confined sample. As is seen from Eq. (|25[) . the energy- 
dependence of the LDOS at fixed position contains ex- 
actly the same information as the position-dependence at 
fixed energy. Therefore, the independence of n(x; Ep) 
at fixed x does not contain more information than the 
density profile at fixed Ep, and does tell nothing about 
the QHE. 

The stepwise increase of the electron density seen in 
Fig. [3] results, of course, from the smooth linear increase 
of the model potential, which is unrealistic, since it ne- 
glects screening effects, and is energetically unfavorable 
in the presence of strong magnetic fields, as has been 
emphasized by Chklovskii et alr^ and following work. 



B. Parabolic confinement potential 

Another simple but instructive model, that also does 
not describe screening effects but allows to consider 
closed, laterally confined equilibrium systems and to cal- 
culate the energy eigenvalues and -functions analytically, 
is the model of a parabolic confinement potential, which 
we write as V(x) = (m*/2) tt 2 x 2 . Combined with the 
parabolic potential describing the effect of the magnetic 
field, this leads to the effective potential 



TO 

~2~ 



X) 



^ 2 {x-X) 2 + ^ 2 X 2 ], (32) 



with uj = yJuj 2 -\-Q 2 and X — (lo c /uj) 2 X. Energy eigen- 
values and -functions are immediately read off, and can 
be written as 



E n {X) 



u r n 2 /X\ 2 



(t)' 



25 3 

n 2 tx^- 



where 



H/(r 



Vn,x{x) 



1 Si ( \ 

2 + 2^ 2 \j) 
(uj c /ui)£ 2 , and 



E n (X), (33) 



(34) 



Since £ < £, the parabolic confinement leads, apart from 
a shift of their center coordinates, to a reduced width of 
the wavefunctions. A sketch of the model potential, the 
energy bands, and the squared energy eigenfunctions is 
given in the lower part of Fig. |4j The dotted lines indi- 
cate the extent of the eigenfunctions, which is constant 
within each Landau band and increases with the quan- 
tum number n of the Landau band. The upper part of 




FIG. 4: (color online) Lower part: potential V(x) = 
(m* /2)Q 2 x 2 for Q/u) c = 0.1, resulting energy bands E n (X) 
for n = 0, 1, 2., and squared and shifted eigenfunctions 
at X = 0. Energies are in units of the cyclotron energy 
E cyc = FiLUc and lengths in units of the magnetic length £. 
Upper part: effective filling factor v c a(x) (solid lines) and 
current density j y (x) (dashed lines, jo = g a eui c / (2ti£)) for 
three values of the Fermi energy, see text. 



Fig. 2] shows, for UbT = 0.01huj c and three values of the 
electrochemical potential /i* = Ep, the electron and the 
current density, calculated according to Eqs. (TT2"]) and 
CEO) with f n<x = f([E n (X) - n*}/k B T). On the scale 
of the figure, the shown results cannot be distinguished 
from those calculated for T = and the same Fermi en- 
ergies. 

Transforming the X-integral into an integral over X 
introduces a pre- factor (uj/ui c ) 2 , which increases the den- 
sity of effective center coordinates X and thereby of 
cigenstates in each Landau band. Referring the effec- 
tive filling factor to this enhanced density of states, 
v c fi(x) = 2n£ 2 (uj c /u)) 2 n(x), leads to the results shown 
in Fig. 21 with plateau values equal to integer multiples 
of the spin-degeneracy g s in regions where the Fermi en- 
ergy is well between two adjacent Landau bands. Note 
that in these regions the usually defined filling factor 
v{x) = (u) /u) c ) 2 v c s (x) > g s (n occ + 1) is larger than the 
number of fully occupied Landau bands with n < n occ . 

The plateaus of the density profile are separated by 
broadened steps, which reflect the structure of the wave- 
function of the highest, partly occupied band. Clearly 
the width of the density profile increases with increasing 
Ep, and, due the symmetry of the considered potential, 
the profiles are even functions of position. The current 
density j y (x), which is an odd function of position, is, 
similar to the density, the sum of the contributions of 
all (partly) occupied bands, and is determined by the 
current efl 2 X/u! c carried by the state \n, X), and by its 
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occupation. Of course, the total current in the considered 
thermal equilibrium state vanishes. 

The calculation of the LDOS from Eq. (|15|) is also 
straightforward, with £ = x/l and i£(E) = X±(E)/£ 
we hnd 



D(E; 



JL 
2nt 



x - 

2 Q2 



where 6 n {E) = 9(2E 



ul(£-g{E)) 
Hu\&(E)\ 



i(E), (35) 



f 1)) and X±{E)/l 
Of course D{E; —x) 



±(oj c /n)y/2E/(f£}) - (2n + 1). 

a;) holds. The LDOS is shown for three different 
positions x by the black lines in Fig. [5j Since for in- 
creasing I a; I the contribution to the LDOS come from 




FIG. 5: (color online) LDOS D(E;x) for the parabolic con- 
finement model with Q/uj c — 0.1 at three positions, x/l = 
5, 10, 15. The black lines are without collision broadening, 
7 = 0, Eq. (|35p. The colored lines are with collision broad- 
ening according to Eqs. (J2TJ and for 7 = V / E cS = 0.015 
and 0.03. Energy in units of E e s = hw, position in units of 
4s = 1 Do = m*/{-Kh 2 ), and g s = 2. 

wavefunctions with increasing \X\, i.e., increasing val- 
ues of energy dispersion \dE n (X)/dX\, the contributions 
of the individual bands become broader and the gaps 
between these contributions become smaller, as can al- 
ready be seen from Fig. 2) The gap between the contri- 
butions due to the lowest bands vanishes for x/l > 20, 
i.e. if i\dE n {X)/dX\/m > 0.2 holds for the relevant 
X- values. This condition is similar to that found in the 
linear-potential model. 

Due to the vanishing energy dispersion in the center, 
[dE n (X)/dX](0) = 0, the contributions of all bands to 
D{E; x) become very narrow and (5-function-like for small 
values of x. With increasing distance from the center, the 
contributions become wider and clearly reflect the struc- 
ture of the corresponding wavefunctions, notably their 



zeroes. Apparently in the energy range shown in Fig. [6] 
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FIG. 6: (color online) As in Fig. [5l but for smaller x values 
and smaller damping. The lower three plots are for x — £, the 
upper three for x = 21. The gap regions are omitted. 

all the zeroes of the wavefunctions arc reflected in the 
LDOS for x = 21, but not yet for x = I 

Figures [S] and [5] also show the effect of collision broad- 
ening on D(E;x), which is most important for small |x| 
values, where it completely washes out the internal struc- 
ture of the individual contributions. For larger values of 
| a: | the collision broadening levels off the maxima and 
smears out the zeroes of the individual structures, and 
at large values, when these structures become broad, the 
collision broadening becomes relatively unimportant. 



IV. MODEL FOR INCOMPRESSIBLE STRIPES 

In the screening theory of the IQHE in narrow Hall 
barsi£~— incompressible stripes (ISs) play an important 
role, which separate neighboring compressible regions, in 
which adjacent Landau levels are pinned to the Fermi 
energy, since these ISs offer the possibility of dissipation- 
less current flow through an otherwise dissipative Hall 
bar. To understand the width and the separation of the 
QH plateaus, i.e. of the B-intervals in which the resis- 
tance quantization occurs, it is important to understand 
the conditions under which an IS can carry a dissipation- 
free current. The calculation s 17 ' 18 were based on a local 
model for the conductivity tensor <7i oc (a:), which was ob- 
tained from the density-dependent conductivity tensor 
<r(n e i) of a homogeneous 2DES of density n e \ by replac- 
ing this density by the local density n(x) of the inho- 
mogeneous system, <t\ oc (x) = a(n(x)). On ISs of finite 
width with constant integer filling factor the components 
of this a\ oc (x) have the quantized values, and dissipa- 
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tionless transport is obtained, if the current is restricted 
to these ISs. It has been argue d 18 ' 19 that the width of 
an IS must be sufficiently large, e.g., more than several 
times the spatial extent of typical wavefunctions near the 
edges of the IS, since otherwise wavefunctions from oppo- 
site sides of the IS would overlap and lead to quasi-elastic 
scattering across the IS. Such QUILLS processes^ would 
lead to dissipation, so that too narrow stripes between 
neighboring compressible regions cannot support the re- 
sistance quantization. 

A suitable quantity containing quantitative informa- 
tion about the ability of an IS to carry dissipationless 
current should be the LDOS in the IS. To calculate this 
LDOS, we have to model the IS with some care. Even if 
the potential within the IS might be well approximated 
by a linear position dependence, at the interesting ener- 
gies around the Fermi energy (i.e. the electrochemical 
potential) there exist nearby states of the compressible 
regions, which may contribute to the LDOS when its gap 
near the Fermi energy becomes small. 
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FIG. 7: (color online) Upper part: filling factor v(x) 
and current density j y (x), in units of jo = (g s /2ir) eui c /t, 
lower part: corresponding potential V(x), energy eigenval- 
ues E n (X) and squared eigenfunctions at center coordinates 
X s /£ = 0, ±2.46, ±4.92, shifted upwards by E n (X s )/huj c . 

To get an idea how the LDOS changes with the width 
of an IS, we consider as a crude model the sum V(x) = 
Vs(x) + Vr(x) of a smoothened step potential 



1. 



Vs(x) = hw c 



^e+Co) 

1 2 



0, 



i < -Co, 

-fo < £ < -f+ , 
-£+■ < £ < £+, 

i+ < t < £ , 



(36) 



and a weak linear potential Vr(x) = 0.1fcsT£, with £ = 
x/l. We take £o = (1 — 7)/ a ) £+ — l/ a - an d k = a 2 /(2 — 
47), with < 7 < 1/2, so that the fraction 2j of the step 



height hu! c is bridged by the linear part of Vs (x) and the 
total width of the stripe is 2£o- In the following we take 
7 = 0.4 if we want to avoid sharp kinks in the potential, 
or 7 = 0.5, if we want to avoid the arbitrarily introduced 
smoothening by parabolic potential regions. 

The weak linear term Vr(x) is added to allow for a 
variation of the electron density in the regions \x/£\ > 
£o- The idea is to simulate the situation in compressible 
regions, where self-consistent screening leads to pinning 
of Landau levels to the Fermi energy, accompanied by a 
variation of the effective potential V(x) of the order of 
k B T. 



A. Thermal equilibrium 

First we consider the system without imposed current 
in thermal equilibrium with constant electrochemical po- 
tential /i* = Ep = 1.5fud c , so that for X <C — I the low- 
est energy band Eq(X) approaches Ep from below, and 
for X 3> £ the second band Ei{X) approaches Ep from 
above. Numerical solution of the eigenvalue problem for 
this potential V(x) with a = 0.1 yields the energy bands 
and the electron and current densities presented in Fig. [7] 
The width (~ Ija = 10£) of the IS is large enough to 
allow for an inner region with constant electron and cur- 
rent densities, similar to the gap regions of Fig. [3J which 
shows results for a linear potential with the same slope 
dV/dx = eE x = —0.1huj c /£. The corresponding LDOS is 




FIG. 8: Energy dependence of the LDOS for the five indicated 
positions x. The broadened step potential is given by Eq. (|36p 
with a = 0.1 (and 7 = 0.4, see text). 

sketched in Fig. [8] for five characteristic values of the po- 
sition x. In the center of the IS, at x = 0, one finds the 
same LDOS as for the linear-potential model with the 
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same electric field, see the lower panel of Fig. Q] As the 
position x moves towards the edges of the IS and leaves 
the regime of linear potential, the individual contribu- 
tions of the different bands become asymmetric and nar- 
rower, since the magnitudes of the slopes \dE n {X)/dX\ 
become smaller. 

For positive a;-values the curvatures of potential and 
energy bands become positive, and the low-energy parts 
of the individual contributions are enhanced, while the 
high-energy parts are reduced, just as we found for the 
parabolic confinement potential in Fig. [SJ For negative 
x-values the curvatures become negative and the asym- 
metry of the individual contributions is inverted, with 
reduced low-energy and enhanced high-energy parts. For 




FIG. 9: (color online) Lower part: Potential V{x) = Vs(x) + 
Vt{x) and corresponding energy bands E n (X) for a param- 
eters as indicated (7 = 0.4). For a — 0.1 shifted squares 
of energy eigenftmctions are indicated. Upper part: corre- 
sponding filling factors and current densities, Ef = 1.5-E cy c, 
ksT — 0.01i? C yc, -Ecyc = de- 



positions close to the high-energy edge of the IS (x/£ rj 
—4.9) there are many nearby states with lower energy, 
but nearly no states with slightly higher energies. As 
a consequence, the individual band-contributions to the 
LDOS for such x- values show a sharp high-energy cutoff. 
Similarly, near the low-energy edge of the IS (x/£ ~ 4.9) 
we find low-energy cutoffs. For positions outside the IS 
> £0 = 6) the individual band-contributions to 
the LDOS are extremely narrow and <5-function-like, very 
similar to the bare Landau DOS. 

Since for the screening theory of the IQHE the exis- 
tence of incompressible stripes of finite width is crucial, 
we present in Fig. [5] for potential steps of different steep- 
ness the resulting energy bands and the electron density 
and current density profiles. Density plateaus with inte- 
ger filling factor v(x) — g s arc obtained for a < 0.1. For 
a > 0.2 the potential increase is so steep that it does not 



lead to a stripe of constant filling factor. Then it makes 
no longer sense to address the step region between the 
flat parts of the potential as incompressible stripe. If we 
assume that the critical steepness is close to a = 0.15, ISs 
do not exist if the width of the potential steps is not larger 
than ~ 7£, which is a little larger than the extent of the 
low-energy wavefunctions. Figure [TU] demonstrates how 




FIG. 10: (color online) LDOS at three positions for several 
steepness parameters a. 

the LDOS behaves in this limit. The positions \x/£\ = 1.2 
are for a = 0.2 well inside the stripe in the region of weak 
curvature, and one sees a similar behavior of the LDOS as 
for \x/£\ = 2.46 in Fig. [51 For a = 0.4, on the other hand, 
\x/£\ — 1.2 is at the edge of the stripe and the situation 
similar to that for \x/£\ = 4.92 in Fig. [8j The situation 
in the center of the stripe, at x = 0, is not so easy to 
interpret. For a = 0.2 the contributions of the individ- 
ual bands to the LDOS are broader than for a = 0.1, as 
indicated in Fig. [lOl but they are still separated by well 
developed gaps, although according to Fig.[S]no IS exists. 
If we increase a further, the gaps shrink, but the LDOS 
vanishes at the energies E = Jho c (n + 1/2), even if the 
potential step becomes very narrow. The reason is sim- 
ple: slightly below E = 1.5Huj c there are nearby states 
in the band Eq(X) and slightly above there are nearby 
states of the band E\{X). But as the energy approaches 
E = 1.5haj c , the center coordinates of these states move 
away from X — and the value of their wavefunctions 
at x — becomes exponentially small. 

This behavior of the LDOS in the center of the IS 
makes the definition of an overlap of the contribution 
of adjacent Landau levels as a criterion for the vanishing 
of the gap in the thermal equilibrium situation useless. 
Things change, however, if we consider a situation with 
imposed current. 
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B. Imposed Hall current 

We now consider an externally imposed current along 
the system and assume (nearly) perfect screening. The 
current is accompanied by a Hall potential Vjj (x) , which 
as a consequence of screening is constant in the compress- 
ible regions and therefore must drop over the region of the 
potential step. Self-consistent screening calculations^ 
show that the width of the incompressible stripes is 
changed by the applied current. It becomes larger, if 
applied and intrinsic current have the same direction, 
and the width becomes smaller, if applied and intrin- 
sic currents have opposite directions. Here we will not 
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FIG. 11: (color online) Lower part: Total potential, cor- 
responding energy bands, electrochemical potential, and 
squared and shifted energy eigenfunctions for X s /£ — 
0, ±1.98, ±3.90. Upper part: corresponding filling factors 
and current densities. Parameters a = 0.13, vh = Vjj/fiw c = 
0.2, see text. 

consider these details and make the simplifying assump- 
tion that the total potential Vt t(x) = V(x) + Vh(x) 
and the electrochemical potential H*{x) are constant in 
the compressible regions \x/t\ > £o and vary linearly in 
the stripe region, i.e., we put in Eq. (|36| 7 = 0.5 and 
£0 = £+ = l/(2a0i and thus suppress the quadratic re- 
gion. To be specific, we take Vjj(x) = Vg F(xjl\ £0) with 

f h 

F(Z,to) = \ -a£, |e|<£o, (37) 

as Hall potential, V tot (x) = V H (x) + hw c [l/2 + F(x/£;^o)} 
as total potential, and (i*(x) = 1.5Hu) c +Vh(x) as electro- 
chemical potential. Numerical results for energy bands, 
wavefunctions, electron and current density are presented 
in Fig. ifora = 0.13 and Vg = 0.2haj c . Although 
the potential has sharp kinks near x = ±3.9£, the en- 
ergy bands E n (X) are smooth and the curvatures near 
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FIG. 12: (color online) LDOS D(E; x) for a steep potential 
step, a = 0.3, with and without externally imposed current, 
for five positions x as indicated. Strength and direction of the 
imposed currents are characterized by the corresponding Hall 
potentials vh = V§//kj c with values 0, ±0.1 as indicated. 



X = ±3.9<? become smaller with increasing n. In the 
linear regime near x = the centers of the wavefunc- 
tions <p n ,x(x) are shifted from X to larger values, as ex- 
pected from Eq. (|2"3"|) . Near x = 3.9£, where the poten- 
tial kink can be considered as limit of a positive curva- 
ture, the wavefunctions are somewhat narrower than near 
x = 0. This is immediately understood from Eq. Q34[) and 
a parabolic approximation of the potential near the kink. 
Similarly, near x = — 3.9£, where the potential kink cor- 
responds to a negative curvature, the wavefunction are 
somewhat wider than near x = 0. 

The most important consequence of the imposed cur- 
rent is that the corresponding Hall potential leads to an 
energetic overlap of the high-energy edge of the lowest 
energy band Eq(X) and the low-energy edge of the next 
band E\{X). Thus the situation near x = is similar 
to that in the linear-potential model, and the individ- 
ual band-contributions to the LDOS will overlap, if the 
region of the potential step will become to narrow. If 
the external current is applied in the opposite direction 
to that of the intrinsic current, the sign of Vh{x) will 
change and, instead of an overlap of the bands, a finite 
energy gap between the bands will result. This will lead 
to energy gaps in the LDOS at x = 0, which will remain 
even if the width of the potential step will become small. 

These results for the LDOS are illustrated in Fig. IT2| 
where we consider a steep potential step, a = 0.3, which 
does not allow for a IS with constant density in the step 
region. Without imposed current, vh = 0, the LDOS at 
x = has no gaps, but has very small values in the middle 
En.n+i = [E n (0)+E n+ i(0)]/2 between the band energies 
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FIG. 13: (color online) LDOS D(E;x = 0) for three value of 
the steepness, a = 0.13, 0.2, 0.3 and three values of the Hall 
voltage, vh = Vg/ftuj c = 0, ±1, as indicated. 

E n (X) at X = 0, D(E n>n+1 ;Q) w 0. If a current in the di- 
rection of the intrinsic current is imposed, vh =0.1, these 
values increase drastically and a considerable overlap of 
the contribution due to adjacent bands is observed. If the 
current is imposed in the opposite direction, vh = —0.1, 
well developed gaps occur around the energies E n>n +i, in 
which D(E; 0) vanishes. On the other hand, the density 
profile, which is not shown, changes only very little due 
to the applied current and remains without any flat part 
in the step region for both directions of the imposed cur- 
rent, similar to the profiles shown in the upper part of 
Fig.i] for a > 0.2. 

Figure [TBI demonstrates how gaps and overlap, respec- 
tively, of the LDOS D(E\ x) in the center of the potential 
step, x — 0, depend on the steepness of the potential and 
on the Hall potential. For Vjj < we observe well de- 
veloped gaps, even if the potential step is so steep, that 
no IS exists ( see Fig. For Vjj > 0, i.e. imposed and 
intrinsic current in the same direction, we have a situa- 
tion as in the linear-potential model, and we can consider 
the overlap as a function of the potential steepness, as in 
Fig. [2] Generalizing Eq. ([26j) by 

O n ,n+l = Dn{E n ,n+l', 0)-D n +l (-EW.n+l i 0)/-Dq, (38) 

where the D n (E; x) are the individual band contributions 
to D(E;x) = Y,n D n(E;x) defined in Eq. JI5J). Results 
are shown in Fig. Q3J 



V. REMARKS AND CONCLUSION 

We have simplified and extended previous calculations 
of the LDOS of a Landau quantized 2DES in the presence 
of a constant, unidirectional in-plane electric field ) 23 i 24 



FIG. 14: (color online) Overlaps Ooi and O12 as functions 
of the potential steepness a for different value of the reduced 
Hall voltage Wj, = V^/Hujc, as indicated by the symbols. The 
lines between the symbols are guides to the eye. The heavy 
lines indicate the corresponding results for the linear-potential 
model, shown in Fig. for r\ = a. 

and we have shown that the LDOS in principle is a use- 
ful concept for further calculations, if the occupation 
probability of energy eigenstates depends only on their 
energy eigenvalues but not on other conserved quanti- 
ties. We have also considered the case of a homogeneous 
2DES supporting a homogeneous dissipation-free Hall 
current, which can be described by the standard meth- 
ods of grand-canonical equilibrium^ but not in terms of 
the LDOS. For this linear-potential model we have also 
quantified the overlap of adjacent band contributions to 
the LDOS, which leads to a closing of gaps in its energy- 
dependence and indicates the onset of quasi-elastic inter- 
Landau-level scattering (QUILLS). In realistic situations 
QUILLS will lead to a breakdown of the IQHE, e.g., 
under high externally imposed currents. To get more 
than an indication of the onset of such breakdown ef- 
fects, one should, however, explicitly consider electron- 
impurity scattering under these conditions, which goes 
far beyond the mere calculation of the LDOS of the 
idealized 2DES. Here we have considered only a simple 
phenomenological treatment of collision broadening and 
mentioned that this is not sufficient under strong elec- 
tric fields. Future work on the necessary generalization 
of the treatment of collision broadening, even within the 
frame of the self-consistent Born approximation, seems 
desirable. 

To get a better understanding of the behavior of the 
LDOS in other situations than the simple linear-potential 
model, we have considered a parabolic potential as a 
model for a laterally confined 2DES, which also allows 
analytical calculation of the LDOS, Eq. (|35|) . Whereas 
D(E;x) at x = exhibits one-over-square-root singu- 
larities of the type 1 / \/E — E n (0) at the energies E > 
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E n (0), for \x\ > there are no singularities but, due 
to the increasing electric field strength, the individual 
band contributions to the LDOS become asymmetric and 
broader with increasing \x\. 

Finally we have calculated the LDOS for incompress- 
ible stripes, which are an essential ingredient of the 
screening theory of the IQHEJ£~— We use a simplified 
model of such an IS, which describes the compressible 
regions (CRs) next to the stripe by nearly constant po- 
tentials, and the stripe region in between by a more or 
less linear potential. One expects that, without collision 
broadening, the LDOS in the CRs approaches the singu- 
lar Landau DOS, whereas in the linear-potential regions 
of incompressible stripes the situation should be simi- 
lar to that in the linear-potential model considered in 
sect. IIII Al In the transition regions between nearly con- 
stant and linear potential the results should be compa- 
rable with those for the parabolic-potential of sect. IIII Bl 
Concerning the energy dependence of the LDOS at char- 
acteristic positions x we find these expectations con- 
firmed. However, the closing of energy gaps and the onset 
of overlap of contributions from adjacent bands, which 
we interpreted in the linear-potential model of sect. IIII Al 
as indication for the breakdown of the IQHE, are now 
more subtle and not so easy to interpret. In the center 
of the stripe region at x = 0, the energy gaps become 
smaller as the distance between the CRs decreases and 
at the energies above and below the Fermi energy the 
centers of the energy eigenfunctions move towards the 
center of the stripe region. But exactly at the Fermi 
energy, which separates the lowest energy bands Eq(X) 
and Ei(X), there are no nearby states and D(Ep; 0) ~ 
even if the CRs come so close that, due to the finite ex- 
tent of the wave functions, no genuine IS with constant 
local filling factor v(x) = g s in an x-interval of finite 
width exists. The situation becomes clearer, if one im- 
poses an external current on the system. This leads to a 
Hall potential in the stripe region, which may increase or 
diminish the intrinsic potential variation across the stripe 
region. If the imposed current has the same direction as 
the intrinsic one, the potential variation increases and 
around Ep the two lowest bands overlap energetically. 



Then the situation is as in the linear-potential model of 
sect. IIII Al and the overlap criterion for the breakdown 
of the IQHE can be applied. If imposed and intrinsic 
currents have opposite directions, the potential variation 
across the strip region decreases and at Ep a gap opens 
between the two lowest bands. Then in D(E; 0) an en- 
ergy gap of finite width around E = Ep remains, even if 
the distance between the CRs becomes so small that, the 
stripe region between them can no longer support a dis- 
sipationless current, i.e. support the IQHE. Of course, in 
a real sample both situations occur simultaneously, since 
the intrinsic currents in the stripe regions of opposite 
sides of the sample have opposite directions. Within the 
screening theory of the IQHE one finds that the width of 
the incompressible stripes is different in both situations. 
If imposed and intrinsic currents have the same direc- 
tion, the stripe is wider than in the opposite case , 17 ^ 19 
however, to determine the widths of the stripes and to 
decide whether they can support the IQHE requires an 
involved self-consistent calculation. 

In summary, the LDOS is an interesting concept, is 
easy to evaluate, and can give some hints on possible 
scattering effects, such as QUILLS, which may lead to 
the breakdown of the IQHE. However, to really under- 
stand the IQHE is much more complicated and requires 
non-trivial calculations. One should include the relevant 
scattering effects, which usually lead to dissipative trans- 
port, and find out, under which conditions they become 
ineffective and lead to the peculiar transport phenomena 
observed in the plateau regime of the IQHE. 
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